The influence of quintessence on the motion of a binary system in cosmology 
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We employ the metric of Schwarzschild space surrounded by quintessential matter to study the 
trajectories of test masses on the motion of a binary system. The results, which are obtained through 
the gradually approximate approach, can be used to search for dark energy via the difference of the 
azimuth angle of the pericenter. The classification of the motion is discussed. 
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I. INTRODUCTION 



Several astrophysical observational data have shown that our universe is undergoing an era of accelerated expansion 
[ll-Q. Therefore, in order to explain this bizarre phenomenon, various models of cosmology have been put forward, 
ranging from the simplest one of a cosmological constant to the scalar field theory of dark energy and modified 
gravitational theory as well other hand, taking dark energy for the presumed cosmological component 

is also used, combined with the Einstein field equations, to deal with some local gravitational issues 0, such as 
the relevant features of black holes. Recently, Kiselev has brought forward new static spherically symmetric exact 
solutions of the Einstein equations, either for a charged or uncharged black hole surrounded by quintessential matter 
or free from it, which satisfy the condition of the additivity and linearity [2]. Otherwise, using a linear post-Newtonian 
approach, Kerr et al. considered the orbital differential equations for test bodies of a binary system in the Kerr-de 
Sitter spacetime and gave the elliptically orbital solution [1^. The solutions for parabolic and hyperbolic orbits can 
be obtained via the formulae of In this paper, we study the trajectories of test masses in a binary system under 
the general metric mentioned in Q . Four discrete theoretical values of the state parameter u) from to — 1 are used to 
solve the orbital equations, respectively. Similar to the classic tests of Einstein's general relativity 12], the formulae 
gotten here can be used to explore whether dark energy exists or not along with the improvement upon observational 
techniques in astronomy. By introducing the effective potential, we also discuss the classification of the motion. The 
metric for a black hole surrounded by quintessence is laid out in Section II. We calculate trajectories of test masses 
in detail in Section III. In Section IV, we display the classification of the motion and a brief conclusion follows. We 
use the metric signature (+, — , — , — ) and make G, H and c equal to unity. 

II. THE METRIC FOR SCHWARZSCHILD SPACE SURROUNDED BY THE QUINTESSENTIAL 

MATTER 

To begin with, Kiselev's work is reviewed briefiy. The interval of a spherically symmetric static gravitational 
field is 

ds^ = e^df - e'^dr'^ - [dO^ + sin^ Odip^) , (1) 

where v and ^ are functions of r. Due to the condition of the additivity and linearity, which implies /i + = 0, the 
energy-momentum tensor of quintessence can be written as 

Tt' = T/ - p„ (2) 



To 



(3c. + 1), 



(3) 
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where pq is the density of the quintessential matter while w is the state parameter. By setting /i = — ln(l + /), we get 

r2/" + 3(l+^^)r/' + (3c^ + l)/ = 0, (4) 

whose solutions are of the forms 

fBH=-'^, (6) 
r 

where A and Tg are the normalization factors. Then, according to the relation 

A should be negative for quintessence because the density of energy is positive. As a resulting, the metric for the 
space surrounded by the quintessential matter can be expressed by 

d^'-{^-^-^-^)dt'-{^~^-^--^ydr^-r' {de- + Bin^ Od^^) , (8) 

where M is the black hole mass and Vg = 2M. Obviously, it can be reduced to the Schwarzschild and Schwarzschild-de 
Sitter metrics by A = and lo = —1, respectively. 

III. TRAJECTORY OF THE TEST MASS AND THE MOTION OF A BINARY SYSTEM 

In relativistic dynamics, the contravariant 4D momentum p^^ is defined as — where m is the rest mass. 

Then the covariant momentum, which is more important in dynamics, can be introduced by = g^i,p'^ ■ When 
particles move in the gravitational field, dynamical conserved quantities are determined by the space-time symmetry 
of the metric field. From the metric ([8]), we may see clearly that t and ip are both cyclic coordinates; therefore, 
conserved quantities of the test mass are the t and ip components of p^, i.e. 

/ 2M \ \ dt 

?'" = "^(^i- — -;3:;TTj^ = co'^st., (9) 

2 2 nd'^ 

P3 — ~rar sm 0— — const.. (10) 
dr 

Due to the symmetry, the test mass moves on the symmetric plane formed by the initial velocity vector (3D) and 
center of force. We regard the direction vertical to the orbital plane as the polar axis. Then the motion of the test 
mass satisfies 

Po and P3 are conserved quantities written in the form 

2M X \ dt 

.^g = L, (13) 

where E and L are constants of integration that represent energy and angular momentum per unit mass, respectively. 
Furthermore, the normalization condition for the 4D velocity, g^^u'^u'^ — 1, also provides the equation 

2M A \fdty / 2M A VV^^^V 2 f d^pV _ . 
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Equations dTTI) to (fT4)) are four first integrals of the geodesic equation, which compose a set of self-contained differential 
equations of particle dynamics. 

After being cleared up, these equations can lead to the relativistic extension of the Binet formula of Newtonian 
mechanics, 

— +u=-+3Mu+ u + ^ u , (15) 

where ^ 7^ 0, because the orbit is not a circle and the second term on the right-hand side is a relativistic correction, 
while the last two are quintessential contributions. We have introduced the dimensionless variable u = GM /r, and 
for convenience u= \/r has been used in the process of the calculation. 

Note that u = GM/r is a small quantity. When w < 0, m also appears in the denominator. So we now evaluate the 
magnitude of the last two terms of Eq. (|15[). First, we make the transformation a = — 3w and /3 = A^''"; then Eq. 



1511 turns to 

" = i^ + 3*^«' + ^(r) +^(-) (16) 



The transformation parameter a satisfies a G [0,3], while uj € [—1,0]. As for /3, we know that the state parameter 
w = —1 represents the standard ACDM model, and therefore j3 ~ A^/'^. From recent cosmological observations, the 
magnitude of the cosmological constant is A w Hq^ /c^ w 10^^-^m^^ with iJp ~ 70kms^^ Mpc"^. As a result, the 
magnitude of ^ is /3 ~ (IQ-^^^i/s ^ iQ-i^. Further more, due to the definition u ~ GM/r, we have the expression 

u GM ^ ' 

By way of a practical evaluation, we take the orbital radius of Mercury, r = 5 x 10^" m, which is the smallest in the 
solar system and for the other parameters we get GM — 1.5 x lO'^ m. As a result, Eq. (IT7)) yields 

B 10^" X 5 X IQ-i" ^„ _ 

^ w 3.3 X 10^1" < 10"^. 18) 

u 1.5 X 10-^ 

Furthermore, from the Binet equation in Newtonian mechanics, 

(Pu M 

whose solution is 

M 

u{ip) = —{l + ecoBif), (20) 

we find that M/V^ ^ u. Therefore, when a > (oj < 0) the last two terms of Eq. ([T5|) are far less than u and Mjl? 
and can be regarded as small corrections. 

According to the linear perturbation scheme put forward in (loj . the last three terms could be treated as pertur- 
bations being a function of kp by substituting for u the unperturbed solution. Due to the supernova's dimming, the 
dark energy should satisfy w < —2/3 to fit the observations [Ts", Ti*]. But here, without essential loss of generality, we 
evaluate the quintessential matter from a theoretical point of view with the state parameter in the range of a; € [—1,0]. 
So we discuss the influence of quintessence on the trajectory of the test mass in four cases. 



A. case of a; = 



When w = 0, Eq. p5| reduces to 



where u is a small quantity, resulting in a slight correction of m^. So we ignore it in the first place, and we find the 
zeroth-order approximate solution 

M + ^ 

Mo = — TT^ (1 + e COS ip) . (22) 
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We presume that the bound orbit is an elhpse with eccentricity e < 1. Via substituting uq into the right-hand side of 
Eq. (|2ip . we can get the first-order approximate solution 



M - 



L2 



1 + e cos tp + 3 [ M 



X\ M + 



eip sm ip 



(23) 



where small higher-order quantities, which are impossible to measure very accurately, are neglected for their insignif- 
icant effects. Furthermore, Eq. (|22|) means that (M -f ■|)/L2 ~ u, and in the weak-field regime, M <C r. Introducing 
e = 3(M -|- ^)[(M -|- ^)/L'^], namely £ <C 1, makes cose(p ~ 1 and sineip ~ e(f. Then Eq. (f23| reduces to 



u{ip) 



L2 



[1 + ecos (1 — e) ip] 



(24) 



When considering precession, we may get the difference of the azimuth angle of the pericenter 



Atp ~ DTT [ — 



/AAf + A2/4 



where we ignore a quantity of e'^, and the second term on the right-hand side belongs to quintessence. 



When w = -1/3, Eq. ^ reduces to 



B. case of oj = —1/3 



— +1-Au= — + MIu^ 
q!(^2 



Similar to the case of w = 0, the zcroth-order approximate solution is 

uo = 1 + e cos {p>^l — 
then the first-order approximate solution reads 



M 



1 + e cos Wl - A + 



3M' 



e(/3 sin ^(ysV 1 — A^ 



Introducing e = makes Eq. (|28l) reduce to 

u (((S) = 1 1 + e cos 
Furthermore, the difference of the azimuth angle is 

M2 



(25) 



Aiy9 = 67r 



L2 (1 - A) 



3/2 



27r 



l-\/l~A 



(26) 



(27) 



(28) 



(29) 



(30) 



Here, is ignored, and quintessence also affects the first term, because the zeroth-order approximate solution is 
changed for the component ip. 



C. case of oj = —2/3 



When w = -2/3, Eq. ^ reduces to 



M 

l2 



A 

22.2^2 ■ 



(31) 
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Obviously, the form turns out to be complicated, because v? appears in the denominator. But as mentioned before, 
A ~ A, so the last two terms are both perturbation ones. Here we could deal with it by the linear perturbation scheme 
in where the motion of test bodies in the Kerr-de Sitter spacetime was studied. Substituting the zeroth-order 
approximate solution uo = | (1 + ecos(^) with p = + -1) for u into the r.h.s. of Eq. ([3T|). we can divide the 
resulting linear equation into two parts: 



--^ + Ml = MIp^^ (1 + e cos C3) . 



(32) 



(PU2 , _i 
— + M2 = p 



dif'^ ' ~ ' 2L2 (1 + ecos(p)' 
After being transformed, Eq. (|33|) can be integrated once to the form 

q C 



d_ r 
di 



+ 



2n3/2' 



(33) 



(34) 



where ^ = cos(p, U = U2 — p ^, q = — §^ and C is a constant of integration. If once more we integrate it, using 
formulae (2.264), (2.266) and (2.269) given in then we get 



1 -e2 

where S" is a constant of integration and 



(1 + e^* sin If) cos (p 

e 



+ C cos if + S sin ip, 



1 / e + cos if 

w = . arcsm 



1 + e cos (f 



(35) 



(36) 



Combining with a particular solution of Eq. (|32|) . the final solution for this case is u — U + p ^ + 3Mp '^eipsinip, 
which can be rewritten in the form 



u{ip) 



(1 + e^" sin if) cos 



M 

-^ + - ) {1 + ecos [(1-e) ^]} 



(37) 



with e = 3M + -1) introduced. Here u may reduce to the result of Schwarzschild case [la] as A — > 0. 



D. case of u 



When w = — 1, Eq. ([T5|) reduces to 



dip 



M 2 A 



(38) 



The unperturbed solution is uo = ^ (1 + e cos (/s) with p = 
result (49) in the approximate solution is 

q 



u{ip) 



3 (1 + e^sin^s) - 



1 + e cos ip 



Similar to the case of w = —2/3, according to the 



M 

_{l + ecos[(l-£)^]}, (39) 



2e H — ) cos If 



where q = —Xp^/L^, e = 3M^/L^ and ^ is represented by Eq. ((36|). When A — 0, this result reduces to the 
Schwarzschild solution. It is different from the two former cases in that the difference of the azimuth angle of the 
pericenter hardly comes out when uj — —2/3 or —1 due to the complicated form. 



6 




FIG. 1: The effective potential for timelike geodesies for variable L with A = and M = 1. Solid line for L > 4, dash line 
for L = 4, dash-dot line for L = 2\/i. 



IV. DISCUSSION AND CONCLUSION 

Here, we turn to a discussion of the classification of the motion. First, we introduce the relativistic effective potential 



r r 



Thus, the radial equation derived from can be rewritten as (^) ~ E"^ ~ V^, which implies that the possible 
motions only occur when ^ V^. Expand Eq. (|40l) . 

1/^ = 1-^ + ^-^-^-^. (41) 



r 



Because the effect of A on is very small, the curves of are almost the same as that in the Schwarzschild space 
[l6j . So we illustrate the effective potential for variable L with A = in Fig. [T] 

From Fig. [1] we see that for L > 4, there is a potential barrier whose peak value is larger than 1 when r is very 
small, while there is also a potential well when r is large. Thus, the possible motions can be classified as three kinds 

E^ < 1, the bound state; 

1 ^ E^ < V„^jaxi the scattering state; 
E'^ ^ ^axj the absorbing state. 



With decreasing L, the central potential barrier decreases in altitude. In the range of 2^/3 < L < 4, V^ax ^ 1j which 
indicates that the appearance of the scattering state is impossible, with the bound and absorbing states being left 
over. If L < 2-\/3, both the peak and the hollow of the effective potential will disappear, and only the absorbing state 
is left over. But note that the above analysis is based on one assumption that the radius of the gravitation source is 
so small that the exterior of the source is applicable for very small r. 

In summary, we have studied the trajectories of test bodies on the motion of a binary system in the presence of 
quintessence. There are four cases of the state parameter uj of quintessence, and we obtain four corresponding orbital 
equations, where it is assumed that the test mass moves round the center of force with eccentricity e < 1. The 
difference of the azimuth angle of the pericenter can be cast in formulae in the two former cases, while the other two 
complicated ones can hardly be done, considering their different forms. Moreover, the effect caused by dark energy 
is too small to be detectable in the solar system. But the common feature of our four cases is that the resultant 
solutions reduce to the outcomes in Schwarzschild space naturally as the parameter A vanishes. 
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